Possible ferro-spin nematic order in NiGa 2 S4 
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We explore the possibility that the spin-1 triangular lattice magnet NiGa2S4 may have a ferro- 
nematic ground state with no frozen magnetic moment but a uniform quadrupole moment. Such 
a state may be stabilized by biquadratic spin interactions. We describe the physical properties of 
this state and suggest experiments to help verif y this proposal. We also c ontrast this state with a 
'non-collinear' nematic state proposed earlier bv lTsunetsueu and Arikawal for NiGa2S4. 
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I. INTRODUCTION 

Spins on a lattice often develop some periodic order at 
low temperatures. It has been appreciated for some time 
that such ordering may be killed even at zero tempera- 
ture by quantum fluctuations. This is particularly true of 
spin systems in low dimension and/or the presence of ge- 
ometric frustration of the magnetic interactions. Indeed 
one of the earliest suggestion of a "resonating valance 
bond" spin liquid state was for the nearest neighbor an- 
tiferromagnetic Heisenberg on the triangular lattice. 
Although, the ground state of this system is now be lieved 
to be the "120° ordered" state fsee lCapriotti et aLp l| and 
references therein), the experimental search for low spin 
quantum magnets on a frustrated perfect triangular lat- 
tice has continued. 

Recently, iNakatsuii et al.p j reported an experimental 
realization of an insulating spin-1 quantum magnet on 
a perfect triangular lattice in the chalcogcnidc NiGa2S4. 
This is a layered material with alternating Ni-S and Ga- 
S planes. Each Ni 2+ ion is in a (distorted) octahedral 
environment of six S 2 ~ ions, and the Ni layer has Ni 2+ 
ions arranged in a perfect triangular lattice. The system 
is an insulator; the only important electronic degree of 
freedom being the spin on the Ni 2+ ion which has an 
electronic configuration of t^e^ (due to the crystal field 
splitting brought about by the octahedral environment 
of the S 2_ ions), and due to the Hund coupling in the e g 
sector, has 5=1. The compound shows unusual mag- 
netic properties. The high temperature susceptibility is 
Curie- Weiss like (with a Wiess temperature of —80 K), 
consistent with an S = 1 moment on the Ni 2+ ions. How- 
ever at around the Weiss temperature, the susceptibility 
begins to drop smoothly, and reaches a finite value at 
T = indicating the absence of a spin gap. The magnetic 
specific heat shows two humps, one at around 10 Kelvin, 
and another broad hump at about 80 Kelvin (the Weiss 
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temperature). However there is no sign of any singular 
or discontinuous behavior suggesting a phase transition. 
The magnetic entropy per spin shows a plateau between 
15 to 50 Kelvin at about a third of the high temperature 
entropy suggesting large degeneracy of low lying excita- 
tions above the ground state. The specific heat, up to 10 
K, shows a T 2 power law suggesting the presence of lin- 
ear dispersing gapless modes in the Ni planes. Moreover, 
this behavior of the specific heat was found to be robust 
to non-magnetic Zn substation in the compound^. In- 
triguingly despite this no periodic spin order is detected 
in powder neutron scattering, and a broad (possibly in- 
commensurate) peak is seen with a wavelength that is 
roughly twice that of the well known 120° spin order 
on the antiferromagnetic triangular lattice. Furthermore 
there is no indication of any lattice distortion associated 
with possible development of spin-Pcicrls order on cool- 
ing to low temperature. 

Motivated by these experiments, 

iTsunetsueu and Arikawajll | have recently proposed 
a spin nematic ground state for NiGa2S4. In such a state 
there is no average magnetic moment ((S) = 0) but 
there is a non-zero quadrupole moment characterized by 
a non-zero average of 
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where S" is the a-spin component operator at site i, S a/3 
is the Kroneker delta symbol. In the state proposed in 
Q , the quadrupole order parameter has a three sublattice 
structure. Precisely if one writes 
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then the "director" n points along three orthogonal di- 
rections in the three sublattices of the triangular lattice. 

In this paper we propose an alternate spin nematic 
state as a possible ground state of NiGa2S4. The partic- 
ular state we propose may be dubbed a 'ferro-nematic' 
in the sense that the average quadrupole moment is spa- 
tially uniform in the ground state. The director vector 
n is independent of site i. This is thus distinct from the 
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Tsunetsugu-Arikawa (TA) statepj which may be dubbed 
a non-collinear nematic. We discuss a number of different 
properties of the ferro-nematic state that makes it attrac- 
tive as an explanation of the properties of NiGa2S4. We 
also point out differences with the TA state that may be 
used to distinguish them in experiments. 

Microscopically the simplest nearest neighbor antifer- 
romagnctic Hamiltonian for a spin-1 triangular magnet 
is expected to have a 120° spin order which is not seen 
in NiGa2S4. Thus the Hamiltonian must apparently in- 
volve other terms that destabilize the magnetic order. 
This is further supported by the observation that bulk 
MS2 is close to a Mott transition and hence may have 
significant charge fluctuations even when insulating. In 
NiGa2S4 the magnetic layers consist of N1S2 sheets - an 
effective spin-only description of these may then plau- 
sibly have sizable interactions beyond the simplest near 
neighbor antiferromagnetic exchange. 

Motivated by these considerations we will discuss the 
possibility of spin nematic order in the framework of the 
model Hamiltonian 

H = H a + H a (3) 

H = jY^Si-Sj-KYjiSi-Sj) 2 (4) 

(id) (id) 

Here Hq is the part of the Hamiltonian that is isotropic 
in spin space. H a refers to small anisotropy terms which 
may be important in pinning any possible long range 
order. Si is the spin-1 operator at the site i of a 
two dimensional triangular lattice, and J, K > 0. Bi- 
quadratic interactions of some strength K will in general 
be present in spin-1 magnets and are known to favor ne- 
matic ordering. 0, In the present case where the bilin- 
ear exchange J is frustrated the effects of the biquadratic 
exchange are somewhat enhanced. We study the ground 
state of such a Hamiltonian in the mean field approxima- 
tion and construct a zero temperature phase diagram. 
We find that at values larger that a critical value of the 
bi-quadratic exchange K/J ~ 1.15, there is a first or- 
der transition from a magnetic(120° order) state to a 
ferro-nematic ordered state. It is the latter that we pro- 
pose to be ground state of the NiGa2S4. We then study 
the expected properties of the ferro-nematic by including 
fluctuations beyond the mean field. 
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where 5 runs over all the neighbors of the site i. 



For the ferro-nematic state the mean spin at any site 



TA also invoked biquadratic interactions to stabilize 
their non-collinear nematic state. However, they chose 
K < 0. Microscopically a biquadratic term may arise 
from a few different sources. First in an underlying two- 
band Hubbard type description of the e g orbitals, the 
spin Hamiltonian arises describes the gain of kinetic en- 
ergy due to virtual fluctuations in a Mott insulating state. 
In this large-C/ perturbation theory the leading term is 
just the usual antiferromagnetic exchange. At higher or- 
ders a variety of terms will be induced which include the 
biquadratic exchange. The sign is readily guessed by just 
considering two sites. In a total spin-2 configuration of 
the two sites, there can be no net gain of kinetic en- 
ergy - so K must be such as to make total spin-2 a high 
energy state. This requires positive K. In some systems 
biquadratic interactions may also be induced due to spin- 
phonon coupling. Indeed if an optical phonon couples to 
the spin bilinear in a bond, its effects can be approx- 
imately accounted for by integrating it out in favor of 
a biquadratic term. Again this leads to a positive K. 
We note that positive K tends to favor the ferro-nematic 
state. 



II. MEAN FIELD AND FLUCTUATIONS 

We now proceed to obtain an approximate ground 
state of the isotropic part of the Hamiltonian Hq adopt- 
ing a mean field approach. Using the definition the 
Hamiltonian Hq can be rewritten as 

H = (J + ^)J2 S i- Si Kj2QfQf, (5) 

(«) («) 

ignoring constant terms. The K term is thereby seen to 
favor ferro-nematic order. For small K we expect that 
the stable ground state will have spiral magnetic order 
which may give way to the ferro-nematic state as K is 
increased. We therefore calculate the mean field energies 
of spiral magnetic and ferro-nematic states to obtain the 
mean field phase diagram. 

We introduce the expectation values of the spin and 
nematic operators respectively as (Si) and (Qf ) at each 
site and obtain the local Hamiltonian at site i as 
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is zero, and the (Q a ^)i has the form 
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with n independent of the site index i. The mean field 
value of q is readily calculated to be —1 corresponding to 
an energy per spin —2K. This energy must be compared 
with that of the spiral state. 

For the spiral state the mean spin at the site i is taken 
to be of the form 



(Si) = msi 



(8) 



with Si a unit vector defined at the site i with position 
vector r, as 



s l = cos (q ■ n) e x + sin (q ■ n) e y 



(9) 



where m and q are, as yet undetermined, magnitude and 
wave vector respectively, and es are two orthogonal basis 
vectors in the plane of spiral ordering. The spin ordering 
will induce a non-zero (Q,f) of the for m 
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with Q satisfying 



1 < Q < !• The four variational 



parameters m,q x ,q y ,Q are determined by minimizing the 
ground state energy of the mean field Hamiltonian JBJ . 

Fig-Ola) shows the variation of the ground state energy 
as a function of K/ J. The standard spiral state with 120° 
three sublattice order is obtained as the ground state 
for (K/J) c < §f ~ 1.143. Beyond this value the fcrro- 
nematic state is found to have lower energy. A study 
of the spin order parameter m (fig. ^b)) shows that it 
jumps at the transition point indicating a first order tran- 
sition. 

The structure of the fluctuations beyond the mean 
field in the spiral state arc well-known and consist of 
three gapless spin wave modes. The ferro-nematic state 
also breaks global spin SU(2) symmetry despite the ab- 
sence of a frozen magnetic moment. Indeed the director 
n picks out a single direction in spin space. Thus this 
state will have gapless Goldstone modes corresponding 
to slow transverse fluctuations of the director. A straigh- 
forward calculation shows the existence of two such gap- 
less modes (corcsponding to the two independent trans- 
verse directions in which the director can tilt) centered 
at wavevector k = 0. At low energies the dispersion is 
linear lo = c\k\ with velocity 
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III. 



PHYSICAL PROPERTIES OF THE 
FERRO-NEMATIC STATE 



We now consider various physical properties of the 
ferro-nematic state with an eye toward interpreting ex- 
periments in NiGa 2 S4. In thinking about experiments 
it is important to allow for the possibility of small spin 
anisotropics that pin the nematic order parameter. The 
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FIG. 1: (color online) (a) Variation of the mean field ground 
state energy E g per site with K/J. Note the "level crossing" 
at (K/J) c = 1.143. (b) The variation of m with K/J. A first 
order phase transition occurs at (K/J) c = 1.143 at which 
m vanishes, (c) The variation of q with K/J. There is a 
jump in q at (K/J) — + , and another at the critical value of 
(K/J) c = 1.143. 
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most important of these is single ion anisotropy that se- 
lects out a particular direction in spin space. Specifically 
consider 



H a = Dj2(Slf 



(12) 



with D > 0. The natural choice is to have the hard 
axis point perpendicular to the Ni-S layers. Note that 
this anisotropy couples directly to the ferro-nematic or- 
der parameter. Thus a non-zero D even if small will 
generally pin the director to point perpendicular to the 
layers. Then in the ferro-nematic state even though there 
is no magnetic order the spins will predominantly fluctu- 
ate parallel to the layers at low temperature. 

In the presence of a non-zero D the linear dispersion 
of the small \k\ director fluctuation modes will be cut-off 
- instead a gap of order D will develop as k — > 0. For the 
specific heat this implies that 



C ~ T 



(13) 



so long as T > To ~ D. In NiGa2S4, precisely this 
behavior is seen in the range 0.35 K < T < 10 K. Thus 
the ferro-nematic state can explain the low temperature 
specific heat provided D is smaller than about 0.35 K. 

Let us now consider the magnetic susceptibility. In a 
single domain sample assuming that the director is fixed 
the pure ferro-nematic state has strongly anisotropic 
magnetic susceptibility. The system is spin-gapped for 
magnetic fields aligned precisely along the nematic di- 
rector, but has finite susceptibility for fields applied per- 
pendicular to the director. A straightforward calculation 
shows that the susceptibility tensor (per site) is 
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where n a is the direction of the nematic director (now 
ferro-aligned at all sites), g is the gyromagnetic ratio, 
and [Ib is the Bohr magneton. If the director is not fixed 
however and is free to rotate it would prefer to orient it- 
self in a direction perpendicular to the external field, and 
the susceptibility would be a constant. In the real sys- 
tem the detailed behavior therefore depends sensitively 
on the competition between the pinning due to the spin 
anisotropy D and the reorientation energy in a field. At 
very low fields the pinning will win and the susceptibility 
will be anisotropic. At higher fields the pinning will be 
overcome and the susceptibility will be a constant inde- 
pendent of the field orientation. As the pinning energy 
per site ~ D while the reorientation energy ~ B 2 the 

crossover occurs at a field strength Bo ^ J — ~ v 7 KD. 

V X 

Therefore the energy scale associated with the depinning 
field Bo is much larger than D since K 3> D {K is ex- 
pected to be of the order of 100 K and D < 0.35 K), and 
thus Bo expected to be a measurably large field. 

In a single crystal sample, the susceptibility at a 
measuring field B <C Bo will therefore be strongly 



anisotropic while at higher fields a constant susceptibil- 
ity independent of field orientation will be obtained. In 
polycrystalline samples (where the existing experiments 
have been done), due to the existence of domains with 
differing orientations of the pining axis, the susceptibil- 
ity will be independent of field orientation even at low 
fields (< Bo)- However as the field increases beyond Bo 
domains with pinning direction parallel to the field will 
reorient their directors thereby increasing the net magne- 
tization. We therefore expect significant increase in the 
magnetization on field strengths on the scale of Bo- 

The constant low-T susceptibility and the T 2 specific 
heat of the ferro-nematic are consistent with the experi- 
ments. However these properties are also shared by the 
TA state. What experiments may distinguish the two 
states? A useful and direct signature is the polarization 
of the spin fluctuation spectrum. In the ferro-nematic 
state (with director perpendicular to the Ni-S planes) the 
spins should primarily fluctuate in the directions along 
the planes at low temperature. In the TA state on the 
other hand there should be no such strong preference for 
the spins to lie along the planes. Thus spin polarized 
neutron scattering experiments on single crystals may be 
able to determine which (if either) of these two spin ne- 
matic states is realized in NiGa2S4. Another useful ex- 
periment that could distinguish the ferro-nematic state 
from the TA state is the direction dependent susceptibil- 
ity measurement in single crystals. As noted above a van- 
ishing susceptibility is expected for a field smaller than 
Bo along the pinning direction in the ferro-nematic state, 
with a significant increase in the susceptibility when the 
field exceeds Bo- Field in a direction perpendicular to 
the pinning direction would obtain a constant suscepti- 
bility independent of the field strength over field scales 
of order Bo and above. On the other hand, in the TA 
state the anisotropy in the susceptibility is expected to 
be much smaller for small fields, and the susceptibility is 
expected to remain constant with increasing field. 



IV. CONCLUSION 

In this paper we have proposed that NiGa 2 S4 may 
have a ferro-nematic ground state characterized by a uni- 
form quadrupole moment. We discussed some of the 
experimental properties of this state and suggested ex- 
pe riments to distinguish it from the alternate proposal 
of iTsunetsueu and Arikawal fll] . In a simple model for a 
spin-1 triangular magnet with bilinear and biquadratic 
terms, the latter promotes nematic order. We suggest 
that the naturally expected sign for the biquadratic term 
prefers ferro-nematic order and not the non-collinear ne- 
matic. We note that while many properties of cither 
spin nematic state (specific heat, susceptibility) resem- 
ble that seen in experiments neither proposal directly 
addresses the short ranged incommensurate spin fluctu- 
ations inferred from powder neutron data. This feature 
of the spin fluctuation spectrum presumably depends on 
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the details of the microscopic spin Hamiltonian which are 
not known at present. 

Jus t prior t o submission of this paper 
cond-mat/0605234 Q appeared which also studies 
the ferro-nematic state in the J-K spin-1 triangular 
magnet, and suggests this as a possible ground state in 
NiGa 2 S 4 . 
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